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ABSTRACT 

Recent observations of giant flares from soft gamma-ray repeaters have exhibited multiple ~ 25- 150 Hz 
oscillations. Frequencies in this range are expected for toroidal shear waves in a neutron star (NS) crust, lending 
support to Duncan's proposal that such modes may be excited in these events. This motivates a reassessment 
of how these waves reflect the NS structure and what role the magnetic field plays in setting their frequencies. 
We calculate the eigenfrequencies and eigenfunctions of toroidal oscillations for a realistic NS crust, including 
a vertical magnetic field at magnetar strengths (B ~ 10 14 - 10 15 G). The lowest radial-order mode has a red- 
shifted frequency of as 28 Hz[l(l+ l)/6] 1 ' 2 , with the prefactor depending on the NS's mass and radius, and 
its crust's depth and composition. This mode is independent of the magnetic field for B <Ax 10 15 G, a limit 
much greater than the inferred dipole magnetic fields for these objects. Though this is a good fit to the observed 
oscillations, only rather loose constraints can be made for the NSs' properties because all that can be fit is this 
prefactor (a single parameter). Modes with shorter radial wavelengths are more sensitive to the magnetic field 
starting at B ~ 2 x 10 14 G and have higher frequencies (~ 600-2000 Hz). The discovery of these modes, 
coupled with the oscillations observed thus far, would provide a powerful probe to the NS crustal structure. 
Subject headings: gamma rays: bursts — stars: magnetic fields — stars: neutron — stars: oscillations 



1. INTRODUCTION 

Soft gamma-ray repeaters (SGRs) are identified by their 
short (~ 0.1 s), recurrent bursts of soft gamma-rays (see re- 
view by Woods & Thom pson! 12004ft . Typical bursts reach 
luminosities of ~ 10 41 ergs s" 1 , considerably above the Ed- 
dington limit. All four of the known SGRs have spin pe- 
riods of w 5-8 s, and three have large period deriva- 
tives (with respect to typical radio pulsars) of ~ 10~ 10 s s" 1 , 
implying dipole magnetic fields of B m (5 - 8) x 10 14 G 
(Woods & Thompson 2004). SGRs h ave been successfully 
explained by the "magnetar" mo del {Duncan & Thompson 
ll992tlThompson & DuncanTT9 95) in which these objects are 
isolated neutron stars (NSs) powered by the decay of ultra- 
strong magnetic fields. 

In addition to the shorter bursts, there have been three "gi- 
ant flares" that released ~ 10 44 - 10 46 ergs of e nergy, which 
took p lace 1979 March 5 from SGR 0526-66 dMazets et al. 
f!979 |) 1998 August 27 from SGR 1900+ 14 (IHurlev et al. 
119991 iFeroci et aJH999). and 2004 December 2 7 from SGR 
1806-20 ( IHurlev et al 120051 iPalmer et all2005l) . In addition 
to the large modulations seen at each object's known spin, all 
three of the giant flares showed signs of periodicit ies a t much 
higher frequencies. From SGR 1 806 - 20. Esrael et al.l ( 2005 ) 
reported a strong oscillati on at 92.5 Hz with wea ker oscilla- 
tions at 18 and 30 Hz. Strohmaver & Watts (2005) discovered 
a 84 Hz oscillation in SGR 1900+14 with weaker oscilla- 
tions at 53.5 and 155.1 Hz and a possible additional feature 
at 28 Hz. Finally, there was some hint of a ~ 43. 5 Hz signal 
in the giant flare from SGR 0526-66 llBarat et al.ll983l) . The 
thrilling implication is that these oscillations may be shear 
modes, excited from crustal deformations during the flare 
dDuncanl|1998l) . If true, they would be the first modes ever 
seen from a NS, initiating a new era where seismology is used 
to learn about NS crusts. 

By extrapolating the calculations of McD ermott et al.l 
(1988) with better estimates for a NS mass, radius, and shear 
modulus, Duncan (1998) provided a frequency estimate for 
toroidal shear modes of a magnetar. Magnetic effects were in- 



cluded by assuming that it acts as a tension adding to the shear 
modulus isotropically. Strohmaver & Watts (2005) showed 
that this estimate matches SGR 1900+ 14's frequencies for 
1 = 2,4,7, and 13 (where I is the angular quantum number). In 
this present work we calculate the modes' dependence on the 
properties of the NS and its crust, and we provide an analytic 
frequency estimate (eq. (9))- We investigate how and when a 
vertical magnetic field plays an important role in determining 
the toroidal frequencies. 

In ^3 we present equations for toroidal oscillations in a 
plane-parallel layer, and in 5j3]we estimate the expected mode 
properties. Numerical calculations of the modes using a real- 
istic crust model are in §4. We conclude in Sj5]with a summary 
our work and a discussion of future studies. 

2. EQUATIONS FOR TOROIDAL SHEAR WAVES 

We assume a NS with mass M = 1 .4M Q and radius R = 
12 km. The top of th e NS's crust is at the depth where 
T = Z 2 e 2 /ak B T w 173 (iFarouki & HamaguchilU993l) . where 

Z is the charge per ion and a = (3/4™,)'^ is the average ion 
spacing, with «, the ion number density. This typically occurs 
at a density 

*,-23*«f t a**(h)\ff(g), (1) 

where T% = T/10 8 K, A is the number of nucleons per ion, 
and we use Fe as the fiducial nucleus. This depth is sim- 
ilar to that found for accreting NSs, rather than like iso- 
lated ones (p top ~ 10 5 g cm" 3 ), because magnetars are young 
(~ 10 3 yrs) and hot from magnetic field decay and recurrent 
flares (Thompson & Duncan 1996). At p k. 4 x 10" g cm" 3 
neutron drip begins, at which point A increases rapidly with 
depth and neutrons play an increasingly important role in pro- 
viding pressure. Th e crust/core interface is at p rj 10 14 g cm" 3 
with A ~ 300-400 ( iDouchin & Haensel200ll) . below which 
the NS is liquid once again. For this reason, the crust can 
be thought of as a cavity where shear oscillations can pre- 
side. Since the pressure scale height at this depth, H = 
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P/pg ~ 3 x 10 4 cm, is much less than the radius, we use a 
plane-parallel geometry with constant gravitational accelera- 
tion, g = GM/R 2 w 1.3 X 10 14 cm s" 2 . We use a Cartesian 
coordinate system with z as the radial coordinate. 

We focus on shear waves, which are c lassified as toroidal 
oscillations ( Aizenman & Sm everslll977l) . These are incom- 
pressible with no vertical displacement, so we can write a La- 
grangian displacement as £ = £ x x + £ y y, with V • £ = 0. We 
assume that the magnetic field is constant and oriented in the 
vertical direction (B = Bz). Though this is a drastic simplifi- 
cation compared to the tangled fields that may exist within a 
magnetar's crust, this provides a simple first step for combin- 
ing shear and magnetic effects. 

Wit hin the crust, the shear modulus is (Stro hmaver et alJ 
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This provides a shear stress tensor iLandau & Lifshitzl 19701) . 
which for toroidal displacements simplifies to be, 



5a u = fi 



dxj dx. 



(3) 



Assuming an oscillatory time dependence, £ oc e lujr , where w 
is the mode frequency in a frame on the NS surface, momen- 
tum conservation in the x- and y-directions is, 



A 5(7 ,. + I WxS)i , 



(4) 



where we have used the Einstein summation convention for 
repeated indices. The perturbed current is 

,2 

(5) 



V X SB-^ri x B 

c 1 



where 6B = V X X B) is the perturbed magnetic field. 
Combining the above equations, we derive the wave equation 
for magnetic toroidal shear modes, 

d£\ B 2 
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where = d 2 / dx 2 + d 2 / dy 2 , and V 2 L | = -/(/+l)//? 2 | in this 
plane parallel-geometry. 

3. ANALYTIC ESTIMATES 

We now make some analytic estimates for the expected 
mode properties. A dispersion relation can be derived from 
equation (|6j by taking d/dz—> ik z and Vj_ — > ik± to be the 
vertical and transverse wavenumbers, respectively. Assuming 
p, w constant, this results in 



(7) 



= v 2 (k 2 + k 2 ± ) + v 2 A k 2 + 



where v s = (/i/p) 1//2 and va = B/{Airp) l l 2 are the shear and 
Alfven speeds, respectively. For low radial-order modes k\ 3> 
ui 2 /c 2 , so we drop the last term for our estimates. 

Equation has two limits, both of which are also found 
numerically. The first is when k± ^> k z , as expected for 
the lowest radial-order mode (which has no nodes in the 
radial direction so that k z ~ 0). Its resulting dispersion 
relation is ui 2 = v 2 ^. We evaluate this frequency at the 
crust/core interface, using a fiducial composition for isolated 



NSs (iDouchin & Haenselll2001l) . The shear modulus at this 
depth is 

. 2 

^ bot = 1.2 x 10 30 ergs cm'Vjf f — 
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where pu = p/10 14 g cm 3 and X n is the fraction of neu- 
trons. Setting k\ = l{l+\)/R 2 and including gravitational red- 
shifting from the NS surface, the observed frequency is 
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where Mia = M/1.4M Q and Rn = R/12 km, and n denotes 
the number of radial nodes for this mode. This frequency 
estimate is within m 5% of what we find numerically. Notice 
that it is strongly dependent on the crustal composition, but 
weakly dependent on the density. 

The fact that this mode is independent of the magnetic 
field is clearly a result of our simple field geometry, so it 
would be useful to have an estimate for the effects of a field 
at an angle a from vertical. This would interact with the 
transverse wavenumber as the shear modulus does, giving 
up- w (y 2 + v\ sin 2 a) k\. Using equation (|8j, the critical mag- 



netic field such that v s 
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When B > Z? cr ; t , we cannot trust our frequency estimate in 
equation l|9}. 

The other important limit is when k z k± as is expected 
for all higher radial-order modes (n > 0) since they have k z ~ 
1 /H. As long as / < R/H, this results in 



(11) 



which shows that such modes are independent of I and sen- 
sitive to B. The vertical wavenumbers of these modes are 
given by the WKB limit, i.e. J k z dz ~ nir where the integral 
is across the crust. This means that B can have an impor- 
tant effect when va > v s in the outer crust, even if va -C v s 
at the crust base. It also shows that there should be a fre- 
quency ratio between the n = and n > modes on the or- 
der of k x /k± ~ R/H ~ 20-30 (as confirmed by our numeri- 
cal ca lculations, also see lCarroll et alJl9 86: McD ermott et alJ 
1988). The observation and measurement of this ratio could 
therefore be used to constrain the depth of the crust. 

4. NUMERICAL CALCULATIONS 

We nex t perform numerical calcu lations using the crustal 
models of Haensel & Pichon ( 1994) abov e neutr on drip (p ss 
4 x 10 11 g cm' 3 ) and lDouchin & Haensel feOOll) for p > 4 x 
10 11 g cm" 3 . Since observed magnetars are young (~ 10 3 yrs) 
and hot these cold, beta-equilibrium models may not be en- 
tirely accurate. An important future calculation would be to 
understand the magnetar crustal structure as a function of age. 
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FIG. 1 . — The shear speed, v s (solid line), and the Alfven speed, va, using 
B = 3 X 10 14 (dot-dashed line), 6 X 10 14 (dotted line), and 1 15 G (dashed 
line). We use the crustal compositions of Haensel & Pichon 1 1994) for p < 
4x 10" gcnT 3 and lDouchin & Haensel <200ll) for o > 4 X 10 11 gcnT 3 . 



The most relevant quantities for our calculations are v s and va, 
which we plot in Figure [2 The key point is that at the base 
of the crust v s 3> Va, even for a magnetic field considerably 
larger than the dipo le fields of B (5 - 8) x 10 14 G inferred 
for these SGRs (Woods & Thomps onl2004ll . 

We set £/R = 1 and the horizontal "traction" to zero in 
the ocean (at p w 5 x 10 7 g cm" 3 ) and integrate equation (|6) 
deeper, shooting for zero horizontal traction at the crust base 
to find u> (the traction includes both shear and magnet ic con - 
tributions, for a complete discussion see Carroll et al. 1986). 
The top boundary condition is equivalent to setting SB = 0, 
and its placement was found to not change our conclusions. 
The bottom boundary condition implies that the crust is de- 
coupled from the core, a poor approximation if the magnetic 
field is as strong as B cr [ t . Future studies should investigate 
what limit this places. We make sure that the horizontal trac- 
tion is continuous at the ocean/crust and neutron drip bound- 
aries. 

Figure |2 shows the eigenfunctions for the three lowest or- 
der modes with B = 2x 10 12 G (left panels) and 10 15 G (right 
panels), all using 1 = 2. The absolute value of the displace- 
ments are plotted, so cusps correspond to nodes. The n = 
mode has no nodes and a red-shifted frequency of 27.4 Hz, 
very close to equation 0. The n > modes are combina- 
tion Alfven/shear waves. When B < (4-irp) 1 / 2 at the top of the 
crust, the crust acts as a hard boundary, trapping Alfven waves 
in the ocean. As B increases, this boundary becomes irrele- 
vant and the modes go all the way to the crust base. Thus the 
n = 1 trapped ocean Alfven wave at B = 2 x 10 12 G becomes a 
shear wave at B = 10 15 G. We also plot the energy density per 
logarithm pressure, which provides intuition of the depths for 
which the mode is most sensitive, 

— — = 27rR 2 l(l+ Du??P/g. (12) 
dlogP 

The n = mode's energy is concentrated at the crust base, so 
it is most sensitive to this depth, explaining why our analytic 
estimate is so accurate. The n > modes have dE/dlogP 
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FIG. 2. — The eigenfunctions and energy density per logarithm pressure 
(measured in ergs) of the three lowest order modes with 1 = 2. The left panel 
displays £ = 2xlO'~G and modes with red-shifted frequencies of 27.4 (« = 0, 
solid line), 521.5 (n = 1, dashed line), and 637.2 Hz (n = 2, dotted line). The 
right panel displays B = 10 15 G and modes with observed frequencies of 27.4 
(n = 0, solid line), 742.8 (n = 1, dashed line), and 1329.2 Hz (n = 2, dotted 
line). 



more constant across the crust, so they are sensitive to when 
va > v s in the outer regions. 

The energy density allows us to test whether the excitation 
of these modes is plausible. A giant flare releases ~ 10 44 - 
10 45 ergs of energy, while from Figure |2] we estimate that ~ 
10 49 ergs is needed to excite the n = mode to £ ~ R. Since the 
mode energy scales oc £ 2 , 10 43 ergs is enough energy for £ ~ 
0.001/? ~ 10 3 cm, which is of order the pressure scale height 
at the top of the crust. It seems likely that a small amount of 
the energy from a flare can excite appreciable amplitudes. In 
contrast, the n > modes require orders of magnitude more 
energy for a similar size, perhaps explaining why they have 
yet to be detected. 

In Figure we plot the red-shifted frequencies of the n = 
0, 1,2, and 3 modes as a function of B. This confirms that the 
n = mode is not affected by the magnetic field and scales 
oc [/(/ + 1)] 1 ' 2 , while the n > modes are independent of / 
and exhibit lu oc B when they are Alfven-like in character. 
This happens at low B for the trapped ocean Alfven waves 
(as described above) and when B>2x 10 14 G. The frequen- 
cies of the trapped oscillation are likely sensitive to our top 
boundary condition, and not robust (whereas the frequencies 
for B>7xl0 12 G are accurate because these modes are con- 
centrated within the crust). 

The numerically calculated n = 0, I = 2,4,7, and 13 modes 
have frequencies of 27.4,50.1,83.8, and 151.0 Hz, respec- 
tively, between m 0.2-6% of SGR 1900+ 14's frequencies. 
This confirms Strohmayer & Watts' identification, and thus 
is strong evidence for a shear wave interpretation of the gi- 
ant flare oscillations. This does not provide strong constraints 
on physical properties because the observed frequencies scale 
<x[l(l+\]] x l 2 . Therefore, only the prefactor can be fit, provid- 
ing a single measured parameter, and as equation (|9j shows, 
this parameter is degenerate with respect to multiple proper- 
ties of the NS. The n > modes would be useful for get- 
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FIG. 3. — The observed frequencies of the n = (solid lines), 1 (dashed 
line), 2 (dotted line), and 3 (dot-dashed line) modes as a function of B. We 
plot I = 2,4, and 7 for « = 0. The « > modes are independent of / (at least 
at these low /-values). Frequencies with B > fi cr j t (shaded region) cannot be 
trusted if the magnetic field has a non-vertical component. Notice that at low 
B avoided crossings occur between the trapped ocean Alfven waves and shear 
waves. 



ting better constraints on the crustal properties, and Figure 
|3 shows that searches should be conducted in the range of 
~ 600-2000 Hz. Due to the energetics arguments described 
above, the powerful SGR 1806-20 giant flare is most promis- 
ing for such a discovery. The 30 and 92.5 Hz oscillations from 
SGR 1806-20 naturally map to / = 2 and 7, respectively, but 
its 18 Hz oscillation has yet to be identified. Though a crustal 
interface wave natural ly provides a frequency in this range 
dPiro & Bildstenl2005l) . magnetic effects cannot be ignored at 
the top of the crust, so further investigation is required. 



5. CONCLUSIONS AND DISCUSSION 

We have calculated the properties of shear modes with a 
vertical magnetic field and compared them to the giant flare 
oscillations seen from SGRs. The n = mode matches the 
observed frequencies and is largely independent of the mag- 
netic field strength. This is not surprising given that the two 
dominant oscillations from SGR 1900+ 14 and SGR 1806-20 
are only 10% different (84 and 92.5 Hz, respectively), while 
the inferred magnetic fields from these systems differ by a 
factor of 1 .4. Some tuning of the interior field strength is re- 
quired if magnetic effects are the predominant factor setting 
the observed frequencies. It seems more likely that differ- 
ences in these NSs' masses, radii, or crusts (as highlighted in 
eq. 0) are the primary reason for the different frequencies of 
these two systems. Nevertheless, future studies should con- 
sider more realistic field geometrie s since they can act on the 
modes in surprising ways (Messios et al. 2001). 

Even if the magnetic field does not play a primary role in 
setting the n = frequencies, it is likely crucial in allowing 
the oscillations to be observed. Toroidal oscillations are in- 
compressible so they are normally decoupled from the ocean 
and do not produce temperature and intensity variations. The 
magnetic field changes things; it creates trapped ocean Alfven 
modes and causes the shear waves to have considerable am- 
plitudes throughout the ocean. A calculation of th e transmis- 
sion t hrough the ocean and magnetosphere (e.g., Bl aes et all 
1989) would be key for understanding the emission mecha- 
nism. Relating this to the observed coherence of the oscil- 
lations, Q ~ 20 - 50, could provide an additional, powerful 
constraint. 
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